On the Polyakov loop in 2+1 flavor QCD 
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We study the temperature dependence of the renormalized Polyakov loop in 2+1 flavor QCD 
for temperatures T < 210 MeV. We extend previous calculations by the HotQCD collaboration 
using the highly improved staggered quark action and perform a continuum extrapolation of the 
renormalized Polyakov loop. We compare the lattice results with the prediction of non-interacting 
static-light hadron resonance gas, which describes the temperature dependence of the renormalized 
Polyakov loop up to T < 140 MeV but fails above that temperature. Furthermore, we discuss the 
temperature dependence of the light and strange quark condensates. 

PACS numbers: 



I. INTRODUCTION 

At high temperature strongly interacting matter un- 
dergoes a transition to a new state characterized by de- 
confinement and color screening (see e.g. Refs. [l], 0] for 
recent reviews). The Polyakov loop is an order parame- 
ter for deconfmement phase transition in SU(N) gauge 
theories. After proper renormalization it is related to 
the free energy of a static quark Fq [H, More pre- 
cisely, it can be defined through the difference in the 
free energy of a system containing a static quark anti- 
quark (QQ) pair at infinite separation and a system with- 
out static charges at the same temperature Foo(T), i.e. 
L ren {T) = exp(-F 00 (T)/2T) = exp(-F Q /T) Q. In 
the confined phase F^ — oo since the static Q and Q 
cannot be separated to infinite distance. Consequently, 
the Polyakov loop which transforms non-trivially under 
the center of the gauge group is zero. In the deconfincd 
phase the static quark and anti-quark could be separated 
to infinite distance due to color screeningwhich means 
breaking of the center symmetry Z(N) [1,1a]. Dynamical 
quarks explicitly break the Z(N) symmetry of the par- 
tition function and Foo is finite since the static Q and 
Q can be now separated to infinite distance by creating 
a dynamical quark anti-quark pair from the vacuum, a 
phenomenon often called string breaking. In 2+1 flavor 
QCD with the quark masses realized in nature there is 
no phase transition related to deconfinement. Moreover, 
the renormalized Polyakov loop cannot be related to the 
singular part of the free energy density 0, Hj]. However, 
the renormalized Polyakov loop is sensitive to the color 
screening in hot QCD medium, at high temperature it 
is closely related to the Debye screening mass (see e.g. 
Q). In the opposite limit of very low temperatures Fq 
is related to the binding energy of a static-light meson 
(see e.g. [13]). Thus the renormalized Polyakov loop is a 
good probe of the hot strongly interacting medium. 

On the lattice with temporal extent N T the renormal- 
ized Polyakov loop is calculated according to the follow- 
ing formula 



cxp(-cA+/2) ( ,Tr J[ U (x ,x) 



(1) 



where Uq (xq , x) is the gauge link variable in the time 
direction and c is the lattice spacing dependent normal- 
ization constant that ensures that the static potential cal- 
culated on the lattice has a certain value at a chosen dis- 
tance [jjj. In the recent past the renormalized Polyakov 
loop has been calculated on the lattice in 2+1 flavor QCD 
with physical quark masses using improved staggered 
fermion formulation 0, [H, ITlTtlil ]. Furthermore, using 
the stout improved staggered action continuum results 
for the renormalized Polyakov loop have been presented 
[l4| . The aim of this paper is to study the renormalized 
Polyakov loop in the low-temperature and transition re- 
gions and to perform an independent continuum extrapo- 
lation using the highly improved staggered quark (HISQ) 
action [15]. While the temperature dependence of the 
Polyakov loop in QCD at high temperatures is very sim- 
ilar to its temperature dependence in pure gauge theory, 
this is not the case for the low temperature and the tran- 
sition regions. To understand at which temperature color 
screening effects set in, it is important to clarify to what 
extent the temperature dependence of the Polyakov loop 
can be understood in terms of hadrons. As mentioned 
above, at very low temperatures the dominant contribu- 
tion to Fq is given by the lowest static-light state. As 
the temperature increases, more massive states will con- 
tribute as well and also the interactions of static-light 
hadrons with the medium will become more important. 
For the description of the bulk thermodynamic quanti- 
ties it turns out that interactions between hadrons can be 
taken into account by adding the contribution of hadronic 
resonances. It is reasonable to assume that the effects of 
interactions of static-light hadrons with the hadrons in 
the medium can be accounted for by adding excited (res- 
onance) states. Therefore, we calculate the renormalized 
Polyakov loop in the approximation of non-interacting 
gas of static- light hadrons and hadronic resonances as has 
been su gge sted recently [l||. Contrary to Ref. [l6| (see 
also [T7I Il8l] ) , where the experimental spectrum of heavy- 
light (strange) hadrons was used together with different 
model considerations, our analysis is largely based on the 
lattice QCD calculations of th e sp ectrum of static-light 
and static-strange hadrons [l9|, |20( . We also consider dif- 
ferent quark model analyses of the heavy-light (strange) 
hadron spectrum, compare them with each other and the 
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available lattice calculations, and use them to estimate 
the contribution of higher lying excited states to L ren . 

The rest of the paper is organized as follows. In section 
II we present our numerical results for the renormalizcd 
Polyakov loop. In section III we discuss the spectrum of 
static-light hadrons and the calculation of the Polyakov 
loop using the hadron resonance gas approximation for 
static-light hadrons. In the course of our numerical in- 
vestigations we also calculated some quantities which are 
sensitive to the chiral transition in 2+1 flavor QCD and 
the corresponding results are presented in section IV. Fi- 
nally, section V contains our conclusions. 

II. NUMERICAL RESULTS 

The chiral and deconfining aspects of the QCD tran- 
sition have been studied by the HotQCD collaboration 
using lattices with temporal extent N T = 6, 8 and 12 
and combination of the tree-level improved gauge action 
and HISQ action in the quark sector Q. This combina- 
tion of the gauge action and quark action was referred to 
as the HISQ /tree action in Ref. Q but here we refer to 
it as the HISQ action for simplicity. For reliable contin- 
uum extrapolations we need at least three lattice spac- 
ings. Therefore we performed calculations using 40 3 x 10 
lattices, using as in the earlier work the rational hybrid 
Monte-Carlo algorithm [2l|. The algorithmic details of 
dynamical HISQ simulations can be found in Ref. (22j. 
As in Ref. [8\ calculations are performed for the physi- 
cal value of the strange quark mass m s and light quark 
masses m; = m s /20. This light quark mass corresponds 
to the pion mass of 160 MeV in the continuum limit [8J, 
which is slightly above the physical value. However, for 
the Polyakov loop this small difference from the physical 
value plays no role. The parameters of the lattice simula- 
tions including the lattice gauge coupling f3 = 10 /g 2 and 
the strange quark mass in lattice units are shown in Table 
U along with the corresponding temperatures. The last 
column of the table shows the accumulated statistics for 
each P value in terms of molecular dynamics time units. 
The lattice spacing a is determined from the r% parameter 
defined in terms of the zero-temperature static potential 
as 



and we use the value r\ = 0.3106 fm 23]. We use 
the parametrization of the lattice spacing and the quark 
masses as functions of the gauge coupling j3 along the 
lines of constant physics that are given in Ref. 0. The 
P dependent normalization constant c that enters Eq. 
(TTJ) was also taken from Ref. [8J. Since we are interested 
in the low-temperature behavior of the Polyakov loop, 
we also performed additional calculations on 32 3 x 8 lat- 
tices at (3 = 6.125, am s — 0.0966 and f3 = 6.175 and 
am s = 0.0906. These correspond to a temperature of 125 
MeV and 131 MeV respectively. We accumulated 3180 



p 


m s 


T [MeV] 


#TU 


6.341 


0.0740 


123 


7289 


6.390 


0.0694 


129 


4540 


6.460 


0.0642 


138 


4200 


6.515 


0.0604 


146 


4918 


6.550 


0.0582 


151 


4390 


6.575 


0.0564 


155 


3910 


6.608 


0.0542 


160 


4290 


6.664 


0.0514 


168 


5000 


6.700 


0.0496 


174 


4990 


6.740 


0.0476 


181 


4990 


6.770 


0.0460 


186 


4990 


6.800 


0.0448 


192 


5310 


6.840 


0.0430 


199 


4990 


6.880 


0.0412 


207 


4990 



TABLE I: Run parameters for 40 3 x 10 lattices. The last 
column shows the accumulated statistics in terms of molecular 
dynamics trajectories. 



molecular dynamics time units for (3 — 6.125 and 3732 
time units for /3 — 6.175. The numerical results for the 
Polyakov loop are shown in Fig. [TJ To obtain continuum 
results for the renormalized Polyakov loop we first per- 
form a smooth spline interpolation of the numerical data 
for each N T . The errors of the spline interpolation are de- 
termined using the bootstrap method. Then we perform 
a 1/N 2 continuum extrapolations at selected tempera- 
ture values. This is motivated by the fact that the leading 
discretization errors in the staggered fermion formulation 
are proportional to a 2 , and therefore for the renormalized 
Polyakov loop they should scale like (aT) 2 = 1/N 2 . For 
T > 133 MeV we have results for L ren at 3 or 4 lattice 
spacings. Therefore for 135 MeV < T < 210 MeV it 
is straightforward to perform a continuum extrapolation 
at selected temperature values separated by 5 MeV. For 
temperatures below 133 MeV we have only few lattice 
data points for iV r = 8 and 10. Therefore in order to 
perform a continuum extrapolation we assumed that the 
coefficient in front of 1/N 2 is temperature-independent 
and is the same as at T = 135 MeV. With these as- 
sumptions we get very good fits for a set of temperature 
values 123 MeV < T < 132 MeV separated by 3 MeV. 
We also performed linear fits with the coefficients in front 
of 1/N 2 treated as free parameters, i.e. fits with zero de- 
grees of freedom. These give values of the renormalized 
Polyakov loop that agree very well with the ones obtained 
from the previous procedure. The errors on the contin- 
uum extrapolated values of L ren at these temperatures 
are taken to be equal to the largest error of the N T = 10 
results. Our continuum estimates for the renormalized 
Polyakov loop are also shown in Fig. [T] and compared 
with the continuum results obtained using the stout ac- 
tion [3] • The two lattice extrapolated continuum results 
agree with each other, except for T = 140 MeV, where 
our results are larger by two standard deviations. 
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FIG. 1: The renormalized Polyakov loop calculated with the HISQ and stout action. The HISQ results for N T = 6, 8 and 12 are 
from Ref. ||. The continuum stout data are from Ref. [l4| . The filled diamonds correspond to the continuum extrapolation 
for the HISQ action. The right panel shows the closeup of the Polyakov loop in the low-temperature region. 



III. THE HADRON GAS MODEL 

As discussed in section I, at very low temperature the 
free energy of a static quark is largely determined by the 
binding energy of the lowest static-light meson. In addi- 
tion, there are contributions from static-strange mesons 
and baryons with one static quark. Thus, following Ref. 
[l6| , at very low temperature the Polyakov loop is given 
by the contribution of the lowest static-light states 

3L re „ = 4exp(-M/T) + 2 exp(-M°/T) + 

Ei Ei(2/ + l)(2j + 1) exp(-M / B °/T), (3) 

where the spin and iso-spin degeneracies of the meson 
states have been taken into account and the summation 
of all iso-spin (7) states as well as of all light and/or 
strange quark angular momentum states (j) for the low- 
est static-light baryons. Altogether we have contribution 
from 6 mesons and 21 baryons [l6j |. The factor 3 on 
the left hand side of the above equation originates from 
the color normalization in the definition of L ren in Eq. 
([T]) and can be seen from the spectral decomposition of 
the Polyakov loop correlator derived in [24| . Interactions 
with the medium are suppressed at very low tempera- 
tures but start to become more important as the tem- 
perature increases. One may try to include the interac- 
tion by assuming that it can be approximated by reso- 
nances. This assumption seems to work quite well for 
bulk thermodynamic quantities 14, 25 27]. The static- 
light meson contains a divergent self-energy contribution 
which needs to be subtracted. This leaves the mass of 
the ground static meson state undetermined, while the 
masses of all other states are given with respect to the 
mass of the lightest state Ei = Mi — Mo- Therefore we 
can generalize the above equation as follows 

L ren = iexp(-A/T)(4 + 2exp(-£ 7T) + 
En,i,j( 21 + W + !) exp(-*WT)). (4) 



Here Eq is the energy of the lightest static-strange meson 
with respect to the mass of the lowest static-light state 
Mq. The first and the second terms correspond to the 
contribution of the ground state static-light and static- 
strange mesons, while the third term corresponds to the 
contribution of the baryon states and the excited meson 
states. The index n denotes different excited states cor- 
responding to the same values of I and j. The renormal- 
ized Polyakov loop depends on the subtracted mass of the 
lowest static-light meson A. This needs to be adjusted 
to match the lattice data for the Polyakov loop at low 
temperatures, i.e. A should be adjusted to the specific 
scheme used for the normalization of the Polyakov loop 
on the lattice. This matching procedure will be discussed 
in subsection C. In the following two subsections we are 
going to discuss the meson and baryon contributions to 
L ren separately. 



A. Static mesons and their contribution to the 
renormalized Polyakov loop 

Static-light and static-strange mesons are character- 
ized by the angular momentum of the light (strange) 
quark and parity j p . The spectrum of static-light 
mesons consist of approximately degenerate pairs with 
j = \l ± 1/2 1, where I is the orbital angular momen- 
tum. The spectrum of static-light (strange) mesons has 
been studied in 2-flavor lattice QCD by the ETMC col- 
laboration for j up to 7/2 that corresponds to orbital 
angular momentum I = 0, 1,2 and 3 and is denoted by 
S, P± , D± , F± [19( . Also the masses of first static-light 
and static-strange excited meson states for 1/2 - channel 
(first radial or S* state in the ETMC notation ) have been 
calculated [191 ] . To get rid of the divergent self-energy 
contribution the masses of different states are calculated 
with respect to the ground state (S state) mass both in 



4 



the light and the strange quark sectors. These mass dif- 
ferences have been extrapolated to the continuum limit 
and to the physical pion mass and are approximately the 
same for static-light and static-strange mesons, see Table 
5 of Ref. [19|. The errors for the above mass difference 
vary between 12 MeV and 37 MeV. With all the spin and 
iso-spin degeneracies the number of states identified on 
the lattice is 96. To calculate the Polyakov loop accord- 
ing to Eq. (|4]) we need to know Eq, the energy (mass) of 
the lowest static-strange meson with respect to the light- 
est static-light meson. We use phenomcnological consider- 
ations to do so. Consider the spin-averaged mass of the 
ground state charmed (bottom) mesons with strangeness 
S = and S = -1 



3M(D*) + M(D) 



4 




3M(B*) 4 


M(B) 


4 




3M(D* S ) 4 


■ M(D S ) 


4 




3M(B*) -f 


- M{B S ) 



= 1975 MeV, 


(5) 


= 5314 MeV, 


(6) 


- - 2076 MeV, 


(7) 


- = 5404 MeV. 


(8) 



M D 
Mb 
M Ds 
M Bs 



Here we used the values of the charm and bottom meson 
masses from Particle Data Group [28j. We get M o s — 
~M D = 100 MeV and M Bs - ~M~b = 90 MeV. Heavy 
quark effective theory predicts that the masses of heavy- 
light mesons and thus also the above difference should 
scale as the inverse of the heavy quark mass mg. Using 
this and the values of Md and Mb as proxies for the 
charm and bottom quark respectively we get a value of 
84 MeV for the difference of the lowest static-strange and 
static-light meson mass, i.e. Eq = 84 MeV for tuq = oo. 
It is interesting to mention that the value of the strange 
quark mass in M S scheme m s (fi = 2 GeV) = 95(5) MeV) 
is close to the value of Eq . So Eq may be interpreted as 
a constituent strange quark mass. 

To study the contribution of higher excited states we 
will use the D s meson spectrum calculated on the lat- 
tice [29| as well as in a relativistic quark model [30l - l33l |. 
The spectrum of D s mesons has been calculated on the 
lattice using improved Wilson fermion actions [29]. One 
needs to establish a relation between the meson masses 
in the static case and the masses of D s mesons. Heavy 
mesons, D s mesons in particular, are characterized by 
nLj with J being the total angular momentum of the 
meson, L being the orbital momentum and n being the 
radial quantum number. Obviously the j = states in 
the static limit are identified with spin-averaged S'-state 
D s mesons. For finite heavy quark mass the L — 1/2 
state in the static limit splits into two states Ll-i and 
Ll, while the L + l/2 state splits into Ll and Ll+i state, 
e.g. P- becomes IP® and lPi, and P+ becomes lPi and 
IP2. The corresponding splittings, however, are small. 
Such a degeneracy pattern is indeed observed in the ex- 
perimentally established positive parity (P-wave) D and 
D s mesons. Furthermore, the mass difference of various 
D a meson states calculated on the lattice and the spin- 



averaged lowest S-state is in reasonable agreement with 
the mass difference in the static limit discussed above. 
Therefore it is justified to use the D s meson masses cal- 
culated on the lattice as proxies for static-strange mesons. 
Since the above difference is approximately the same for 
strange and light quark cases we can use the same mass 
difference also for the light quarks. That allows to in- 
clude the following excited states into the analysis: 2P, 
2D and 2P. We identify the mass of 2L_ meson in the 
static case with the lowest 2Ll D s meson mass and the 
2L + with the higher 2Ll D s meson. With all the spin- 
iso-spin degeneracies this gives 90 states. 

To include even higher excited states we use quark 
model predictions. The quark model can predict cer- 
tain qualitative features of the heavy-light and static- 
light meson spectrum correctly [30l - l32l |. However, the 
quark models also have problems. In the static limit the 
mass of the P+ state is smaller than the mass of the P_ 
state just the opposite to what is observed on the lattice 
fl9| . The mass of the P + in a quark model, calculated in 
Ref. [32|, is 230 MeV below the lattice result. Similarly, 
the mass of the 2S state is 441 MeV below the lattice 
result [l|| . Comparing the results of Ref. [13, HH to Ref. 
[32! ] one may conclude that the model dependence is small 
for the IP meson states, however, in the B s sector the 
masses of the 25* states differ by about 300 MeV. We try 
to account for these problems by assigning a theoretical 
error to the masses of higher excited states. 

We took the results of quark model calculations of D 
and D s mesons [331 ] to estimate the contribution of 35, 
AS, 5S, 3P and 1G states. From comparison of the re- 
sults of different quark models as well as the comparison 
to the lattice results discussed above we estimate the un- 
certainty of the masses of the higher nS states (n > 3) 
to be 300 MeV, while for the other states we estimate it 
to be 150 MeV. It turns out, however, that contribution 
of these states to L ren is negligible up to temperatures 
of 210 MeV for which the model makes sense. The ex- 
cited states discussed so far should include all the possible 
states up to mass of 2 GeV above the ground state mass. 
It is unlikely that individual resonance states can be ob- 
served above that energy. In Fig. [5] we show the con- 
tribution of meson states to the renormalized Polyakov 
loop. We normalize the results by Lq — 4exp(— A/T)/3. 
The contribution of all static meson states calculated on 
the lattice in Ref. [l9| is shown as the band, the dashed 
line includes the contribution of the higher excited states, 
while the thin solid line corresponds to the first static- 
strange meson only. The uncertainty band is determined 
by the errors of the static meson masses calculated on the 
lattice [19( . We see that at low temperatures the only me- 
son state that has a significant contribution to L ren in 
addition to the lowest state is the lightest static-strange 
meson, though the contribution of excited states cannot 
be completely neglected. Excited states (up to IP) be- 
come very important at higher temperatures, T > 140 
MeV. Finally, the contribution of higher excited states is 
quite small and is only visible for T > 170 MeV. 
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FIG. 2: Meson contribution to L ren normalized by La = 
4 exp(— A /T) /3. The band shows the contribution of all static 
meson states calculated on the lattice (see text). The dashed 
line includes the contribution of higher excited states. The 
thin solid line corresponds to the contribution of the lightest 
static-strange meson only. 



B. Baryon contribution to the renormalized 
Polyakov loop 

The spectrum of baryons with one static quark has 
been studied by the ETMC collaboration [20] in 2 flavor 
QCD. The masses of static baryons with positive and 
negative parity and angular momentum of light quarks 
j = and 1 have been calculated. These states cor- 
respond to the ground state and first orbital excitation 
of baryons with one heavy quark, i.e. to 1/2+ , 3/2 + 
and 1/2 - and 3/2~. Counting the iso-spin and angular 
momentum degeneracies the lowest positive and nega- 
tive parity baryons correspond to 79 states. The cal- 
culations have been carried out at one lattice spacing. 
The lack of continuum extrapolation is not of great con- 
cern, since based on the studies of the static meson spec- 
trum, cutoff effects are expected to be small compared to 
the statistical errors. The results presented in Ref. [2(| 
depend somewhat how the lattice spacing is set. More 
precisely, using /„■ the lattice spacing was determined to 
be 0.079(3) fm, while using the nucleon mass the lat- 
tice spacing turned out to be 0.089(5) fm. In our anal- 
ysis we use the values of the masses obtained by fixing 
the lattice spacing through the nucleon mass TOj\r since 
this procedure gives a value of the ro parameter that 
is consistent with other determinations I8l. Il3l . namely 
r = 0.473 ± Q.09(stat.) ± 0.16(syst.)Jm (34|. Setting 
the scale with / w gives ro = 0.42 fm [35| which is much 
smaller than any other determination. The lowest lying 
positive and negative parity baryons give a fairly large 
contribution to the renormalized Polyakov loop, in fact, 
the largest contribution next to the ground state mesons. 
However, it turns out that higher excited states cannot 
be neglected. We can use quark models to estimate the 



contribution of higher lying baryon states to L ren . 

The spectrum of baryons containing one heavy (c or b) 
quark has been studied in relativistic quark model [36[ 
and in relativistic quark-diquark model [37|. The anal- 
ysis of Ref. [Hj], however, was restricted to A& iC and 
Sb !C baryons. We will use the spectrum of excited heavy 
baryons containing b quark as a proxy for the spectrum 
of higher excited baryon states with a static quark. The 
masses of baryons with a static quark are determined 
by the angular momentum of the light quarks j. There- 
fore the heavy baryons form doublets with almost the 
same mass that correspond to the same angular momen- 
tum j of the light quarks and total angular momentum 
J = j ±1/2. In our analysis we consider the mass differ- 
ence of the baryon with the lower angular momentum in 
the doublet and the spin averaged mass of B(B*) mesons. 
These mass differences obtained in a quark model are 
compared to the static baryon spectrum for the lowest 
positive and negative parity states calculated on the lat- 
tice. It turns out that the agreement between the lattice 
results and the model calculations is quite good if the nu- 
cleon mass is used to set the lattice spacing. In fact, the 
lattice results agree with the model calculations within 
the errors. For the Af, and family we also find good 
agreement between the diquark model and Ref. [36j for 
the lowest states of both parity. In our calculations we 
use the spectrum calculated in Ref. [37j which corre- 
sponds to baryon states with angular momentum up to 
J = 11/2 equivalently to j = 5 of the light quarks and up 
to 5 radial excitations. Counting all the spin and iso-spin 
degeneracies these correspond to 984 states. 

As discussed above, different model calculations agree 
with each other for the lowest positive and negative par- 
ity states. Unfortunately, the agreement is not that good 
for the higher excited states. To estimate the sensitiv- 
ity of the Polyakov loop to the model uncertainty of the 
higher excited baryon states we calculated the contribu- 
tion of excited Aq and Eg baryons to L ren including all 
states up to J = 7/2 using the results of Ref. [36| and of 
the diquark model (37j . The contributions of At, to L ren 
is a factor of two larger if one uses the spectrum from 
Ref. [36[ compared to the case where the A& spectrum 
from the diquark model is used. On the other hand, the 
contribution of the baryons is factor of two smaller if 
one uses the results of Ref. ;36j instead of the results of 
the diquark model. Therefore we estimate that contribu- 
tion of the higher excited baryon states is uncertain by 
factor 2.5. This is the largest source of uncertainty in the 
hadron resonance gas model for T > 170 MeV. The con- 
tribution of baryons to the renormalized Polyakov loop is 
shown in Fig. [3l The contribution of the static baryons 
identified on the lattice in Ref. [2(| is shown as the solid 
line and the band. The error band corresponds to the 
uncertainty that has been evaluated using the errors on 
the static baryon masses. The contribution of all baryon 
states to L ren is shown as the dashed black line. At tem- 
peratures T < 120 MeV the contribution of the baryons 
is below 10%. It becomes significant above that temper- 
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FIG. 3: The contribution of baryons to L ren normalized by 
Lo — 4 exp( — A/T)/3. The solid line and the band correspond 
to the contribution of the lowest positive and negative parity 
baryons, while the dashed line corresponds to the contribution 
of all baryon states. 



ature. The contribution of higher excited states becomes 
significant only for T > 140 MeV. Therefore, as it will be- 
come clear in the next subsection, the uncertainty in the 
comparison to the lattice data due to the excited states 
is small. 



C. Comparison with lattice results 

Let us compare the hadron resonance gas model results 
with the lattice data discussed in section II. The compar- 
ison is easiest in terms of the free energy of an isolated 
static quark Fq(T) — — Th\L ren {T). The renormaliza- 
tion procedure of the Polyakov on the lattice introduces 
a scheme dependence. Therefore, for the comparison of 
the hadron resonance gas with the lattice data one needs 
to adjust the parameter A in Eq. (|4]). Fitting the pa- 
rameter A to all available continuum lattice data in the 
temperature interval 123 MeV < T < 135 MeV we get 
A = 580 ±5 MeV with x 2 /d.o.f = 0.49. The comparison 
of the lattice data with the hadron resonance gas model 
is shown in Fig. |H The solid line and the band cor- 
respond to the hadron resonance gas result with all the 
states discussed above and its uncertainty. The dashed 
line corresponds to the contribution of the ground states 
only. The figure shows that the contribution of excited 
states is significant already for the lowest temperature 
available in the lattice calculations. The hadron reso- 
nance gas model can describe the lattice results on the 
renormalized Polyakov loop up to temperature 140 MeV, 
however, clearly fails above that temperature. All the ex- 
cited states included in the analysis are not sufficient to 
explain the rapid decrease of the static quark free en- 
ergy. This result agrees qualitatively with the findings 
of Ref. [Hi if the stout data for T < 140 MeV are used 



FIG. 4: The free energy of a static quark Fq(T) calculated 
on the lattice and compared with the resonance gas model 
(solid line). The uncertainty of the hadron resonance gas 
model is indicated by the band. The dotted line is the hadron 
resonance gas model result with the ground state meson and 
baryon contribution only. 



to normalize the static free energy in their analysis (c.f. 
Fig. 4 of Ref. [l6|). It is also possible that exotic hadron 
states can explain the discrepancy between the lattice re- 
sults and hadron resonance gas model for T > 140 MeV. 
Another possibility could be the partial restoration of 
the chiral symmetry and the corresponding change in the 
static- light hadron masses, as indicated in a recent lattice 
study [381 ] . 



IV. TEMPERATURE DEPENDENCE OF THE 
QUARK CONDENSATES 

As a byproduct of our study we also calculated the 
quark condensates (i/j^q- The quark condensate needs 
a multiplicative renormalization, and for non-zero quark 
masses also an additive renormalization. It is easy to see 
that the leading additive divergence is proportional to the 
quark mass and is quadratic in the cutoff (inverse lattice 
spacing). Therefore it was proposed to study the follow- 
ing combination, called the subtracted quark condensate 



Ai,.(r) 



(#>i, T -SW>v 

W>J,0-^(^>a,o' 



(9) 



Here q = I and s corresponds to light and strange quarks, 
while the subscript x — 0, r refers to zero and finite 
temperature expectation values. The expectation val- 
ues (4>ip) q , x are normalized per single flavor. Sub-leading 
divergences proportional to the quark mass cubed and 
the logarithm of the cutoff are expected to be small for 
the physical values of the light quark masses. We cal- 
culated A/ iS on N T — 10 lattices. Combining this with 
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FIG. 5: The subtracted chiral condensate calculated with the 
HISQ action. The open diamonds correspond to the N T = 12 
results obtained using the fx scale [8J (see text). 



the published results of the HotQCD collaboration ob- 
tained on iV r = 8 and 12 lattices Q we performed a con- 
tinuum extrapolation. First we interpolated the lattice 
data for each N T using a smooth spline and estimated 
the error of the spline by bootstrap analysis. Next we 
performed a continuum extrapolation at selected values 
of temperature separated by 5 MeV within the interval 
145 MeV < T < 200 MeV assuming a 1/N? behavior. 
We studied the variation of the extrapolated result with 
varying the fit range in N T . These variations have been 
included in our final error estimate. For T < 170 MeV 
the N T — 6 data have not been included in the analy- 
sis as they are incompatible with a 1/iVj? behavior. For 
T < 140 MeV no N T = 12 data are available so the 
extrapolation had to rely on iV r = 8 and N T = 10 data 
only. The numerical results for different N T as well as the 
continuum extrapolations are shown in Fig. [5] The con- 
tinuum extrapolated results are slightly above the con- 
tinuum results obtained with the stout action 14]. This 
difference is expected due to the slight difference in the 
light quark masses used in the two calculations, namely 
mi = m s /20 versus m; = m s /27 in Ref. [l4j. In Fig. 
[5] we also show the iV r = 12 HISQ data obtained using 
the lattice spacing determined from the kaon decay con- 
stant fx @|- These data are systematically above our 
continuum estimate. 

Alternatively, we can get rid of the ultraviolet diver- 
gences in the quark condensate by considering the follow- 
ing combination, which is called the renormalized quark 
condensate I8| 



a r 

g 



d + 2m a r$((ipi>} g , T -{il>Tp) qfi ), q = l,s. (10) 



Here d is a normalization constant that is related to the 
light quark condensate in the chiral limit. More precisely, 
d = 2m s rf('tptp) ifi(mi —> 0). With the values of m s and 
(tpip)l,o(mi -> 0) from Ref. [H we get d = 0.0232244. 



The quantity defined in Eq. ([10]) is closely related to 
the renormalized quark condensate {tpip}}i introduced in 
Ref. [Hj]. Using our N T — 10 results and the published 
HotQCD results for N T = 6, 8 and 12 we perform a 
continuum extrapolation for A^. As for A; :S we first 
perform a smooth spline interpolation and estimate the 
errors of the spline by bootstrap analysis. Then we per- 
form a 1/iV 2 continuum extrapolation for selected values 
of the temperature separated by 5 MeV in the interval 
145MeV < T < 200 MeV based on the interpolation and 
its errors. We performed extrapolations using subsets 
of the available N T values and the differences in the ob- 
tained fit values for were treated as systematic errors 
and entered into our final error estimate. For T < 140 
MeV the continuum extrapolations is based on N T = 8 
and 10 data only. 

The lattice QCD results for Af and A^ are shown 
in Fig. [B] along with the continuum extrapolations. We 
also show the HISQ N T = 8 data obtained using the lat- 
tice spacing from fx in this figure which seem to agree 
quite well with our continuum result, except for T > 180 
MeV, where they are systematically lower. Our contin- 
uum results for A^ are slightly larger than the continuum 
results obtained with the stout action [l4| ■ This is again 
expected to be due to the difference in the light quark 
masses (see discussion in Ref. Q). Finally we would like 
to note the large difference in the temperature depen- 
dence of Ap and Af . The decrease of the renormalized 
strange quark condensate is much more gradual than of 
the light and Af- reaches half of its vacuum value only 
at T ~ 200 MeV. 



V. CONCLUSIONS 

We studied the renormalized Polyakov loop in lat- 
tice QCD using the HISQ action and obtained results 
in the continuum limit for temperatures 120 MeV < T < 
210 MeV. Results obtained with the HISQ action are in 
the a 2 scaling regime for N T > 6. Our continuum results 
agree well with the earlier findings obtained using the 
stout action [l4| . We studied the question of the physics 
origin behind the increase in the Polyakov loop, or equiv- 
alently, the decrease in the free energy of a static quark 
Fq. At sufficiently high temperatures the decrease in 
Fq is associated with onset of color screening which also 
leads to the same decrease in the energy of a static quark 
at leading order [|oj]. For temperatures T < 140 MeV the 
decrease in Fq could be explained in terms of hadron res- 
onance gas model. For larger temperatures the decrease 
in Fq appears to be significantly larger and cannot be 
explained in terms of conventional static-light (strange) 
hadron states. It remains to be seen whether this rapid 
decrease is due to the contribution from exotic static- 
light hadrons or some other mechanism. In the latter 
case its implication for color screening is not clear, espe- 
cially in view of recent lattice results on the static energy 
of QQ pair which do not indicate large significant screen- 
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FIG. 6: The light (left) and strange (right) renormalized quark condensate calculated with the HISQ action. Open diamonds 
correspond to 7V T = 8 HISQ results obtained with the fx scale §. 



ing effects for T < 200 MeV [HI, H|. Currently one of 
the largest uncertainties in the Polyakov loop calculations 
within the hadron resonance gas model comes from the 
excited baryon states. Clearly an improved lattice cal- 
culation of static baryon spectrum would be very helpful 
in this regard. Another open issue is the contribution of 
exotic static-light hadron states. 

We also revisited the calculation of the quark conden- 
sate. Combining our N T = 10 results with published 
HotQCD results we performed continuum extrapolations 
for the subtracted quark condensate as well as for the 
strange and light renormalized quark condensates. We 
found that lattice results for these quantities in the tran- 
sition region follow the expected a 2 scaling if N T > 8. 



numerical simulations have been performed at NERSC 
and on BlueGene/L computers at the New York Center 
for Computational Sciences (NYCCS) at Brookhaven Na- 
tional Laboratory. The smooth spline interpolation and 
bootstrap analysis was performed using the R package. 
We thank S. Mukherjee for his help with the R pack- 
age and F. Karsch for reading the manuscript and useful 
comments. We also thank E. Megias, E. Ruiz Arriola and 
L.L. Salcedo for reading the first version of this paper and 
pointing out the normalization error in the hadron reso- 
nance gas expression for the renormalized Polyakov loop, 
which we subsequently corrected. 
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